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This work presents the derivation of the large time and distance asymptotic behavior of the field- 
field correlation functions of impenetrable one-dimensional anyons at finite temperature. In the 
appropriate limits of the statistics parameter, we recover the well-known results for impenetrable 
bosons and free fermions. In the low-temperature (usually expected to be the "conformal") limit, 
and for all values of the statistics parameter away from the bosonic point, the leading term in the 
correlator does not agree with the prediction of the conformal field theory, and is determined by 
the singularity of the density of the single-particle states at the bottom of the single-particle energy 
spectrum. 
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I. INTRODUCTION 

This is the last article in the series of papers [l], 0, Q in which we study rigorously the large time and distance 
asymptotic behavior of the temperature dependent field-field correlation functions of one-dimensional impenetrable 
anyons. In this work, we present the derivation of the final results for the asymptotics of the time-dependent correlation 
functions. As in the case of "static" (same-time) correlators, for which the asymptotic behavior was computed in 
[H , the starting point of our analysis is the determinant representation for the correlators found in [H, • With the 
help of this representation, we are able to derive a system of differential equations for the correlators, which is the 
same as the one for impenetrable bosons 0,0]) but with different initial conditions. The asymptotic behavior of the 
correlators is computed then by solving the matrix Riemann-Hilbert problem that is associated with the obtained 
system of the differential equations. The most striking feature of the time-dependent asymptotics found in this work 
is the fact that its leading term is non-conformal. It contradicts the predictions of the conformal field theory (or, 
equivalently, bosonization) that were derived for the one-dimensional anyons in @, 0] . This is in contrast to the static 
correlators which do agree with the conformal field theory. 

The model of impenetrable anyons considered in our series of papers is arguably the simplest physical model of 
one-dimensional particles with fractional exchange statistics, and is the anyonic generalization of the impenetrable 
Bose gas first studied by Girardeau Despite its simplicity, the model is closely related to the realistic models 
of transport of anyonic quasiparticles of the fractional quantum Hall effect [9j. The model of impenetrable anyons 
can also be viewed as the infinite-repulsion limit of a more general model of anyons with 5-function interaction of 
finite strength, called the Lieb-Liniger gas of anyons, and suggested in Introduction of the fractional exchange 
statistics in one dimension requires additional convention for the direction of the particle-particle exchanges [7|,l9(. This 
implies that for finite anyon-anyon interaction, the anyonic wavefunction is discontinuous at the coincident particle 
coordinates, the fact that makes the physical interpretation of the finite-interaction case difficult. Nevertheless, the 
Lieb-Liniger gas of anyons can be well-defined mathematically, and has received considerable attention in the last few 



years. As a result of these efforts, we know the Bethe Ansatz solution K | of this model, the low-energy properties and 
the connection with Haldane's [ll[ fractional exclusion statistics [IH, 13 j, the thermodynamics the ground-state 
properties and the low-lying excitations Q. Various techniques were used to study the correlation functions 
(mostly for the physically- motivated impenetrable case) such as: the Fisher- Hartwig conjecture 0, [13) bosonization 
@, conformal field theory fij, numerical calculations 0,(13], and the replica method [Hj]. The present paper, together 
with our previous papers [l|, 0, H[ , is devoted to the exact calculation of the asymptotic behavior of the correlation 



functions using the techniques developed for impenetrable bosons |4|, | a. l3fj. l31l. l32t 13311 . It should be mentioned that 
other models of the one-dimensional fractional exchange statistics |2ll. l22l |23|. l24l |25|. l26l W\ can also be found in the 
literature. In particular, the quantum inverse scattering method with anyonic grading was developed recently in [2S| . 

The main result obtained in this work is the large time and distance asymptotics of the field-field correlator 
(^>{x2, t2)^{x\, t\f)T of impenetrable anyons at finite temperatures. This result can be expressed conveniently in the 
rescaled variables [see Eq. ©] in which 



(^(x 2 ,t 2 )¥(x 1 ,t 1 )} = VTg(x,t,j3,K), 



(1) 
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and the function g(x, t, [3, re) is denned below. To do this, we need to introduce several quantities: 

C(x,t,/3,K) = - J \x~2tX\\-a\ip{X 2 ,P,K)\ dX, I(/3,k) = ( S(j ln<^(A 2 ,/3,re) dXj , (2) 

and 

y(A 2 ,/3,re) = e A^ /3 + 1 » ( 3 ) 

where the branch of the logarithm is chosen so that lim^oo In <,c>(A 2 , /3, re) = 0. With these definitions, our result for 
the function g(x,t, [3, re) can be stated as follows. In the large time and distance limit: x,t — > oo, with x/t — const, 
the asymptotic behavior of the field-field correlation function |T]) is given by 

g(x,t,0,K) = i" 2 /2 e C(x, t ,3,K)+«/(/3,«) [coi -l/2-i Ve 2it(A2+/3) + ^(it^-ixA*) + (1/Vt)] , (4) 

where A s = —x/2t and f = — — In |<p(A 2 , /?, re)|. Other notations are: 

AJ = - (/3 + + ?r 2 re 2 ) V2 /V2 T » + \//? 2 + 7T 2 ^ 2 ) V2 / V2 , (5) 

with re £ [0, 1] being the statistics parameter: k — for bosons and re = 1 for fermions, cq and ci are some undetermined 
amplitudes, and the upper and lower signs correspond, respectively, to the space-like and the time-like regions defined 
by x/2t > \fj3~ and x/2t < \ff3~. Equation ((4|) also assumes the condition \$ty/f3 + inn — x/2t\ > ^s/f3 + inn, where 
one should take the positive branch of the square root, the meaning of which is clarified in the main text. It might be 
argued that for any finite re, the second term in the parenthesis of the asymptotics ((4]) is exponentially small compared 
to the error term and therefore should not appear there. The presence of this term is justified, however, by the fact 
that it becomes dominant in the bosonic limit re — > 0, when SAj — ► 0. It is interesting to note that for all re 0, 
the first, leading term of the asymptotics (|3]) is not the one predicted by the conformal field theory or bosonization 
0, 0], and only the second, sub-leading term gives the conformal part of the asymptotics, as demonstrated explicitly 
in Section IVlTl 

The plan of the paper is as follows. Section [TT] describes the determinant representation for the correlation functions 
obtained in which is used in Section [Ml to obtain differential equations indirectly describing these functions. The 
relevant matrix Riemann-Hilbert problem is introduced in Section [IV! and its asymptotic solutions in the space- like 
and the time-like regions are presented in Sections IVl and IVT1 The complete results for the correlators are summarized 
in Section IVIIl and their analysis in the bosonic, fermionic, and the low-temperature ("conformal") limit is given in 
Sections IVIIII and HXl In the two Appendices we (A) discuss the large time and distance asymptotic behavior of the 
correlators of free fermions, and (B) present detailed analysis of the function C(x,t, (3, re). 

II. DETERMINANT REPRESENTATION FOR THE FIELD-FIELD CORRELATOR 

The second-quantized form of the Hamiltonian of the Licb-Liniger gas of anyons is 

H . / * (ftKCIlMM + =*'M*'M*W*M - . (•) 

where h is the chemical potential and c is the coupling constant, assumed in our case to be infinite to make the anyons 
impenetrable. The anyonic fields satisfy the commutation relations of the usual form 

*{xi)&(x 2 ) = e- l ^ £(xi - x ^^{x 2 )^{x 1 ) + 5(xi - x 2 ) , 



¥(xi)¥(x 2 ) = e inKe{xi - X2) ^(x 2 )&(xi) , 

with e(x) = x/\x\, e(0) = 0. The commutation relations become bosonic for re = 0, and fermionic for re = 1. We are 
interested in the asymptotic behavior of the space, time, and temperature-dependent field-field correlator defined as 

{^{x 2 ,t 2 )^ l {x 1 ,t 1 )) T = Tr e _ H/T • 
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In 0, we have obtained the following representation for the correlator: 

(* (x 2 , h)&{ Xl , ti)) T = (J-G'(t 12 ,x 12 ) + -^j det(l + Vf) 



(7) 



a=0 



where x a {, = x a — xj,, t a j, = t a — tb, a, b = 1, 2 , and det(l + Vr) is the Fredholm determinant of the integral operator 
with the kernel 



Vf(\,iJ.) = cos 2 K2) exp --t 12 (A 2 + + - I12 (A + /i) 



2 y ^ ' 2 
E(X\t 12 ,x 12 ) - E((j,\ti 2 ,xi 2 ) a 



2tt 3 



E(X\ti, 2 ,x 12 )E(fj,\t 12 ,xi 



(8) 



which acts on an arbitrary function /(A) as 

Wf) (A) 



Vf(\,fi)f(p) dp. 



The functions G'(t 12 ,xi 2 ) and E(X\t\ 2 ,xi 2 ) in Eqs. (0 and ([5]) are defined by 



G'(ti 2 ,xi 2 ) 



e 



it\2l*< —1X12^ 



dfi, 



and 



dp + 7rtan(7TK/2)e l * 12A ~ lx " x , 

-co M-A 

where P.V. denotes the Cauchy principal value, and $(A) = $(A, T, h) in Eq. ([8]) is the Fermi distribution function of 
the quasiparticle momentum A at temperature T and chemical potential h: 



1 



l + e (A=-h)/T 



The correlator ([7]) depends on five variables: time, distance, temperature, chemical potential and the statistics 
parameter. It is convenient to rescale three of them and the momentum A by temperature: 

x = (xi - x 2 )VT/2 , t = (t a - h)T/2 , (3 = h/T , X^X/Vt. (9) 

Then the explicit dependence of the correlator on temperature is simple and is given by Eq. ([1}. To see this, one needs 
first to obtain a more manageable expression for the field correlator (J7J). In the rescaled variables ([H]), the functions 
G and E are given by 



G'(t,x) = VTG(t,x), G(t,x) = 



- 2 zt\ 2 - 2 ix\ j, 

e dX , 



and 



E(X\t,x) = P.V. / dfi 



-lit^, — lixti 
jl — X 



+ 7r tan(7TK/2)e" 



- 2 ity -2ix\ 



We introduce the two functions e±(A): 



e_(A) = COs(W2) v/^(A) e * tA2+a * . 



and 



e+(A) - e_(X)E(X) 



cos(7tk/2) 



y^(Xje ltx +lXx 



P.V. 



dfi 



-litti — 2 ixfi 

fj, — X 



7rtan(7r«;/2)e" 



-1it\ A - 2 ix\ 



(10) 



(11) 



(12) 



(13) 
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In terms of these functions, the kernel {8]| of the integral operator appearing in ([7]) is expressed as 

V£(X,») = V T (X,(i)-^A T (X, f i), 

with 

V T (X^) = e+{X)e - {fl) - e - {X)e+M , (14) 
A — /i 

and 

A T (\,fi) = e+(A)e+(/x) . 
In what follows, we also need some basic formulae from the theory of Fredholm determinants: 

°° (-l) n+1 

logdetfl + V) = V - — Tr V n , (1 + VY 1 = 1 - V 1 + V 2 + ... , 

£ — ' n 

n=l 

where y™(A, ^) is determined successively by V n ( A, ^) = J V(X, t^F™ -1 (v, [i)dv and V 1 (A, fi) = V(X,[i). The trace is 
defined naturally as Tr V = J V(X, X)d\. Then, Tr V 2 = J J V(X, n)V(p, X)dXd[i, and so on. Using these relations, 
one can see directly that 

<9det(l + ^; 



= --^Tr [(1 + V T )~ l A T ] det(l + V T ) 



da 

which together with Eq. (|10[) gives Eq. |T]) for the correlator with 

g(x, t, (3, k) = -^e 2ttl3 (tt [(1 + Vt^M - G{t, x)) det(l + V T ) . (15) 

The integral operator Vt whose determinant appears in (fTS")) is of a special type called "integrable" operators 
[J, [H, ■ This type of integral operators have kernels of the "factorizable" structure similar to Eq. (fl4| , and are 
ubiquitous in investigations of correlations functions of integrable quantum systems and distribution of eigenvalues of 
random matrices. If an operator is integrable, the resolvent operator defined as 

R T = (l + V T )~ 1 V T , (1 + V t )(1-Rt) = 1, 
is also of the same type, which means that the resolvent kernel that solves the integral equation 



R T (X,n)+ / V t {X,v)Rt{v^) dv = V T (X,fj,), 

J — oo 

is also factorized as in (fH[) : 

/ + (A)/_(/.)-/_(A)/ + ( At ) 

Kt = ■ 

A — n 

The functions /±(A) are the solutions of the integral equations 



/±(A)+/ V T (X,fi)f ± (ji)dfx = e ± (X). (16) 

J — oo 

Now we can introduce an important class of objects called auxiliary potentials defined as 

r + co 

Bi m (x,t,(3,K) = / ei(X)f m (X)dX, l,m = ±, (17) 



and 

/+oo 
Xei(X)f m (X)dX , l,m = ±. (18) 
-oo 

Due to the symmetry of Vt(X,/j,), we have = B |_. One can see directly that Tr[(l + Vt)~ 1 At] = B ++ . 

Therefore, defining b ++ = B ++ — G, we obtain the following representation for the function (|15l) in the time- and 
temperature-dependent correlator (TTJ): 

g(x, t, /?, k) = --U 2jt/ V+(x, t, 0, k) det(l + V T ) . (19) 

Since g(x,t, (3, k) = g(— x, t, (3, — k) and g(x,t, (3, k) = g*(x, — t, (3, — k), to study the correlator, it is sufficient to 
investigate only the case x > 0, t > 0. 
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III. DIFFERENTIAL EQUATIONS FOR THE CORRELATION FUNCTIONS 

Obtaining the differential equations directly for the correlation functions at finite temperature is an extremely 
difficult task. One can, however, obtain a system of partial differential equations for the auxiliary potentials, and 
show that the derivatives of the logarithm of the Fredholm determinant 

a(x,t,j3,K) =logdet(l + y T ), (20) 

is expressed in terms of a combination of the auxiliary potentials and derivatives. The differential equation for the 
potentials are obtained as follows. First, we define a two-component function 



F(X) 



/+(A) 

f-W J ' 



and look for three matrix operators L(A), M(A), N(A) which depend on the auxiliary potentials and their derivatives 
and satisfy the Lax representation conditions 

L(A)F(A) = , M(A)F(A) = , N(A)F(A) = . 

The differential equations for the potentials are obtained then from the compatibility conditions for the Lax repre- 
sentation 

[L(A), M(A)] = [L(A), N(A)] = [M(A), N(A)] = 

which should be valid for any value of the spectral parameter A. 

Specific calculations follow closely those for the impenetrable bosons [3, @, and their main ingredient are the 
following relations 

d x E(\) = -2iG — 2i\E(X) , 
d t E{\) = -2i\ 2 E -2i\G + d x G, 

d P E{\) = 0, (21) 
d\E{\) = -(4:it\ + 2ix)E-4itG, 

(22) 

which can be proved directly from the definitions (fTU|) and (fTTl) of the functions G and E. Here we only present the 
results of the calculations. 

Lemma III.l. The potentials C(x,t, /?, k) can be expressed in terms of the potentials _B(x,i, (3, k) and G(x,t) as 
follows 

C++ = -^d x B ++ — 2GB + - + — B ++ , 

C— = --d x B— - B+-B— , 
2 



Theorem III.l. Define 



— — C— + — B 2 — B++B- 



mi 



d 



n = g-9+ = b++B — , p = g-d x g + -g+d x g~. 
Then <?_ and g + satisfy the separated nonlinear Schrodinger equation 



id t g+ = 2(ig + + -d 2 x g + + ^g\g- , 



id t g. = 2l3g_ + -d 2 x g- + 4g 2 _g+ , (23) 
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and 

—2idtn = d x p 

The equations containing the (3 -derivatives are 

df3d x g+ _ d p d x g- 
9+ 9- 

and 



= V, 



-id t <p + 4dpp = , d x ip + 8dpn + 2 = 0. 

The previous theorem characterizes completely the potentials B and 6++ (£?++). The other potentials can be 

expressed in terms of these two as 

d x B+_ = 2ib ++ B— , d t B+_ = -p , (24) 
dpB + - = -ix/4 - itp/4 , 

and 

d x (C+_ - C_+) = (B++ - 2G)d x B— - B—d x B++ . (25) 
Finally, we have the following theorem 

Theorem III. 2. The derivatives of the logarithm of the Fredholm determinant a (x,t, j3, k) are given by 
d x a = -2iB+- , 

d t a = -2iGB— ~2i(C+- + C-+) , (26) 
d a = -2itdp(C+- + C-+) - 2ixd l3 B + _ - 2itB—dpB ++ + 2it(B ++ - 2G)9 /3 S__ 
+2(d p B ++ )(d fj B^)-2(d fj B + ^ 2 . 

All the differential equations above do not depend on the statistics parameter, and are the same as those obtained 
for impenetrable bosons in [1, Q . The statistics parameter appears only in the initial conditions which can be extracted 
from the equal-time field correlator studied in @, [2i| . The same phenomenon was noticed also for static correlators 
at T = [17] and finite temperature 0,1111. 



IV. MATRIX RIEMANN-HILBERT PROBLEM 



The discussion in the previous sections implies that with the use of the differential equations, the large-time 
and -distance asymptotic behavior of the field correlator can be extracted from the corresponding behavior of the 
auxiliary potentials. A powerful method of obtaining the asymptotics for the potentials is the formalism of the matrix 
Ricmann-Hilbcrt problem (RHP). Here, we consider a specific matrix RHP associated with the integrable system that 
characterizes the potentials. Solution of this RHP will allow us to obtain the asymptotics of the potentials and field 
correlator. In details, we are interested in finding a 2 x 2 matrix function x(-M) nonsingular for all A 6 C, and analytic 
separately in the upper and lower half planes, which also satisfies the following conditions 

X-(A) = x+(A)G(A), x±(A)= lim x(A±ie), Ael. 

X (oo) = I. (27) 
Here / is the unit 2x2 matrix and G(X) is the conjugation matrix defined only for real A and given in our case by 

G(A) = ( 1 ~ 2 - e +( A )MA) W + (A) \ 

v ' \~2-Kie_(X) 1 + 27T«e + (A)e_(A) J v ' 

The functions e±(A) appearing in this equation are defined in \Y&\ and (|f 2p . The matrix function x(A) depends also 
on x, t, (3, and k, but this dependence is suppressed in our notations. Also, in what follows we will consider k € (0, 1]. 
The case of impenetrable bosons, k = 0, requires a special treatment presented in [3, l33|. 
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A. Connection with the auxiliary potentials 

In this section, we show that in the limit of large A, the auxiliary potentials can be extracted from the solution of 
the RHP (j2"7| . To do this, one can see first (as shown, e.g., in Chap. XV of [5|) that the RHP is equivalent to the 
following system of singular integral equations 

I((11 = r + ifi*ta, ACM. 
2m j_oo /t — X — iv 



Multiplying from the right with 

H(X) 



1 e+(A) 
e_(A) J ' 



and introducing x(X) = x+(A)i?(A), we transform these equations into 

X(A) = H(X) + — / d/i, 

2mJ_ 00 fi-X-iO 

2m J^oo n — X — tu 



where 



and 



G(X) 





2TTie-(fi) 27ri(e+(A)e_ (/i) — e_(A)e+(/i)) J' 



X (X) = x+WHW = ( Xn ' + !xi Xll,+ ffi e+ S t X12A , x! e "S x! ) • (29) 
1 X2i,+ [X) X2i,+ (A)e+(A) +X22,+ (A)e_(A) ' 



The integral equations for X12 and X22 are 

°° Xi 2 (At)(e+(A), 

fi — X 

and 



m , /■ +oc Xi 2 (M)(e + (A)e_( M )- e _(A) e+ ( M )) 
Xi2 = e+(A)+/ d/i, AeR, (30) 



m , /" + °° X22(/i)(e + (A)e_(/i) - e_(A)e + (/i)) 

X22 = e-(A)+/ d/i, Ael. (31) 

J-00 M- A 

Taking into account that the functions /±(A) satisfy the integral equations (fT6|) . where the kernel Vr(A,/i) is given 
by Eq. (TH| , one can see directly from (f3T)]) and (T3Tj) that 

Xi2(A) = /+(A), X2 2 (A)=/_(A). (32) 

Also, Eq. ([29]) gives that xn(A) = Xn,+(A), X2i(A) = X2i,+ (A). Therefore, using Eq. ([32]) we obtain the following 
integral equations for xn,+ (A) and X2i,+ (A) 

J-oc M-A-iO 7- 00 /t-A-iO 

Continuing analytically into the upper half plane, taking the limit of large A, and using the definitions of the auxiliary 
potentials (fTTll and (f!8| . we obtain from these equations 

Xii(A) = l--^- + -^C- + + O^V (33) 
X2i(A) = - ic__ + O (iV (34) 
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In order to obtain similar expansions for Xi2(A) and X22(A), we proceed in the following fashion. First, Eq. (|29[) gives: 

X2 2 (A) = /-(A) = X 2i,+ (A)e + (A) + X22 , + (A)e_(A) . 

Then, using Eq. (|34|) and the large-A expansion of the integral equation (fl6|) defining /-(A) we rewrite this equation 

as 

e_(A) - \ (e+(A)i?__ - e_(A)B_+) - ^(e+(A)C__ - e_(A)C+_) + • • • = 

1 - 1 ~ \ . A . 1 (1) , 1 (2) 



e + (A) ( - -^C — J + e_ (A) I 1 + -x 22 J + ^X 22 J ) + ■ ■ ■ (35) 



Comparison of the two sides of this equation implies that 



X2 2 (A) = 1 + jB + _ - ^C+- + O (± ) (St.) 



The expansion for Xi2(A) can be derived through similar steps: 

Xl2 (\) = l B++ + ±C^ + + o(j-). (37) 

Collecting the results and (pJT)) . we see that in the large-A limit, the auxiliary potentials follow from 

the expansion of the solution of the RHP (|2T|) : 



*«='+K=l:- &) +0 GfJ- 



B. Transformations of the RHP 



It will be useful to perform several transformations on the RHP ([27)1 . The first one is 

x(A)=x(A)xo(A), 

with 

Xo(A)=(^ J j , a(A)=_/_ oo m _ a d/i. 
Using the fact that the boundary values of the function a(A) on the real axis are: 

f+oo -2it[i 2 —2ix(i, 
-2it\ 2 -2ix\ ' 



o±(A) = ±nre~ 2ltA ~ llxX + PV. / dfi , 

J-oo fJi — A 

it can be shown that the matrix x(A) solves the transformed RHP 

X-(A)=x+(A)G(A), Ael; x(oo)=/, (38) 

with (5(A) = Xo+(A)G(A)xoi(A) given explicitly by 

g (x) _ ( 1 - i?(A)(l + e 1 ™) 27ri(i?(A) - l) e - 2rfA2 - 2iAa; 

1 V -~ cos 2 (™/2)tf(A)e 2ltA2 + 2jA * ' 1 - 0(A)(1 + e- 1 ™) 

The specific form of the second transformation depends on whether we are considering the "space-like" (x/2t > ^/J3) 
or the "time- like" (x/2t < \f]3) region. 
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1. Transformation in the space-like case 
As a first step, we need to introduce the functions 

2 „ s e A 2 -/3_ e i7rK 



p(A', /?,«)= eA2 _ /3 + 1 , (39) 



a(A)-exp^--L / +0O ^-l nv ,( M 2 ,/3, K )^ . (40) 



and 

rvf W = pyn <; — 

27T« J-oo M - A 

The latter is the solution of the following scalar Riemann-Hilbert problem (for more information on scalar RHP see, 
e.g., Hi) 

a_ (A) = a + (A)[l-tf(A)(l + e 47rK )], AeR; a(oo) = 1 . 
Then, the second transformation in the space-like case is: 

$(A) = ^(A) e -' T3lnQ ^ , 
where 03 is the third Pauli matrix. The new matrix function $(A) solves the matrix RHP 

*_(A) = $+(A)G*(A), Ael; *(oo)=I, (41) 
with the conjugation matrix G$(A) = e ff3 ha ( A )G(A)e-' r3 lna(A). 

r m _ ( 1 p(A)e- 2 « 2 - 2i ^ \ 

G$(A) " I (? (A)e 2 ^ 2 + 2 ^ 1 + p(X)q(X) ) ' (42) 



where 

— _9-WU mi 2 

=A 2 -/3 



P (A)^-2^[a_(A)] 2 ^ 2 !""^^ (43) 



and 

=,A 2 -t9 



g(A) = -g!cos 2 (W2)K(A)]- 2 . 3 . . (44) 



2. Transformation in the time-like case 

The transformation in the time-like case is similar to the one performed in the space like case. The difference is 
that the function a(A) is now defined as (note the change of the sign of k) 

a(X) = cxp i — - — : f — ^— In </?(/i 2 , (3, —ft) 



2tt« J-oo (i - A 
and is the solution of the scalar Riemann-Hilbert problem 

a_ (A) = a+(A)[l-??(A)(l + e-" rK )], AeR; a(oc) = 1 . 
The new matrix $(A) = x(A)e +CT3 lna ( A ) solves the same RHP (|4"Tj) but now with the conjugation matrix G$(A) 

e -CT 3 lna WQ(X)e +a ' 3 lna(A) . 



where p(A) and g(A) are 



, , _ / l+p(\)q(\) p( A)e- 2 ^ 2 - 2 ^ . 
" \ g (A)e 2! ' A2 + 2,lA 1 1 



A 2 — /3 

p(X) = -2™[a_(A)] 2 ^ — , (46) 



and 

2i I 

- cos 2 (ttk/2) [a+ ( A)] - 2 2 _ 
7r e A p — e 



q{X ) = CO s 2 (7r K /2)K(A)]- 2 -^— — . (47) 
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C. Potentials in terms of the $ matrix 

In Section IIV Al we showed that the auxiliary potentials we can be extracted from the large-A expansion of the 
solution x(A) of the RHP (|'27p. However, since we explicitly will be finding the asymptotic solution of the RHP (I4ip . 
we need to express the potentials in terms of the matrix. The computations necessary to do this are presented 
below only in the space-like case, the time-like case being similar. The first step is to obtain the large-A expansion of 
all the terms in the relation 

*(A) = x(A)xo 1 (A) e -^ lna ^. 

Explicitly, we have in the limit A — > oo: 

1 f -B-+ B++ \ 1 / -C-+ C++ \ . „ ( 1 



^) = /+ a UbI: bZ + * (-cZ. a?- +0 (x* 



y 1 f0 -G\ , 1/0 -G« \ , n ( 1 



™ = I+ A + A + °U 



e -3lna(A) =/+ l f-ao 1 

A \ a I V A 3 / 



where G is given by (fT0|) , and 



1 



+ OC 



a Q = 7T~ I lnip(p 2 ,/3,n) d[j,. (48) 



Considering a similar expansion for ^(A) 

W ~ + A ^ ($1)21 ($1)22 J + \»\ ($2)21 ($2)22 J U 3 
and equating the terms with equal powers of A, one finds 



B+- = -($1)11 -«o, &++ = ($1)12, £— = -($1)21 , 

G + _+G_ + + 5__ G = ($ 2 ) 22 - ($2)11 • (49) 

In the time-like case, similar computations give 

B+- = -($1)11 +"o, &++ = ($1)12 , £—=-($1)21, 

G+_ + G_+ + 5__G = ($ 2 ) 22 - ($ 2 ) u , (50) 

with 

1 / >+ °° 

a = — / lnip(^ 2 ,^, -k) d/x. (51) 

V. ASYMPTOTIC SOLUTION OF THE RHP. SPACE-LIKE CASE 

We are interested in solving the RHP (|41[) in the limit of large x > and t > 0, but with finite ratio x/t — const. 
If one compares the solution with the corresponding solution for the static case (without time t), the analysis in the 
time-dependent case is more complicated due to the presence of the stationary point of the phase 

4>(x, t, A) = tA 2 + xX , 

where d\(f> = 0. This condition gives 

A s = -— ,x>0,i>0,A s <0. 
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The asymptotic analysis of the RHP has to pro perly take into account this stationary point. In this work, we do this 
by employing the method pioneered in l3rl l37j ] , also used for the impenetrable bosons 0, [33[ . The main ingredient 
of this approach is the Manakov ansatz [3y| which provides an approximate solution <I> m (A) to the RHP (j4Tjl . The 
Manakov ansatz in the space-like region is different from the one in the time-like region, even though the results 
obtained from both forms of ansatz will be the same in the leading order. The asymptotic analysis is based on the two 
assumptions: (i) the RHP is solvable, and (ii) the boundary values of < &±(A) on the real axis are uniformly bounded 
in the limit t — ► oo. These assumptions can be proved following the Sections 6 and 7 of (33j . Also, we require that 

\sRy/f3 + inK\ -x/2t > |9f ^(3 + inK\ . (52) 

The meaning of this inequality is discussed below (see Sec. IV C|) . While this condition is not essential in that one 
can analyze other regimes as well, it is always satisfied, in particular, in the more interesting low-temperature case 
(3»1. 

A. Manakov ansatz 

The space-like region is defined by 

\ a <-y/p, = h/T>O, 
condition that can be expressed in more explicit notations as 

(xi - x 2 ) > v F (t 2 - h) > , 

where vf is the velocity of excitations, which in our model of impenetrable anyons coincides with the Fermi velocity 
of free fermions, vf = 2kF, with kp — Vh, in the conventions used in the Hamiltonian (|6|). The Manakov ansatz in 
the space-like region is given by 



where 



P(x) = — f + °° M^-M (M)e - 2l0(a: , t ^ ) (53) 

2tt« J-oo M - A 



± ^ v2. q Qj,yw,w d (54) 
2tt« J^oo M - A 

and the functions p(X) and q(X) defined by Eqs. (|43|) and (|44|) . The function S(\) is the solution of the following 
scalar RHP 

6+{X) = 8^(X)[l+p(X)q(X)ri(X s -X)] ) Ael, <5(oo) = 1 , 
with 77(A) denoting the step function 



??(A) 



1, A > 0, 
0, A < 0. 



This scalar RHP problem can be solved explicitly (see, e.g., [33]), and if we take into account that 

l+p(A)g(A) = |^(A 2 ,/3, K )| 2 , 

the solution is 
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1. Properties of 8(X) 

Before we show that $ m (A) is an approximate solution of the RHP (|4Tj) , it is useful to investigate some of the 
properties of the function 5(X). For A S (A s ,oo), integration by parts gives 

_L f S ^i n |^2 ;/3;K) | 2 = l ln(A _ As)ln | ¥)(A 2 ; ^ K) |_l I 3 ln^- A|d(ln|^( M 2 ,A«)l)rfM- 
2m J_ oa n - A m m J_ oc 

Introducing two quantities: 

v(\ a ,(3, k) = - In \<p(\*,0, k)]- 1 > , (55) 

7T 



1 ^ 



and 

7(A) = - / ln|/i-A|d(ln|^,/3,K)|)d M , 

one can use this relation to write <5±(A) for A € (A s ,oo) as 

5±(A) = (A-A i )^exp(i 7 (A)) J 

where (A — A s )± are the boundary values of the multi-valued function (A — A s ) ll/ defined in the complex plane with 
the branch cut along the ray (— oo, A s ]. When A e (— oo, X s ), integration by parts for singular integrals (see, e.g., [38j . 
pp. 18) gives 

5±(X) = cxp J ± In |^(A 2 , 0, k)\ + - ln(A s - A)± In \cp(\l,0, «)| - -, [ ' In | M - X\d(ln \<p(fj?, 0, «)|)d/x 1 . 

This equation can be rewritten as 

MA) - (A - A.)£ exp(* 7 (A))|^(A 2 , 0, K^V^.fl k)^ 1 . 
in the notations used above. Therefore, the function 5(X) in both regions is 

MA) = (A - \ s )± exp(«7(A)) (|^(A 2 , / 3, K )||^(A 2 ,/3, K )r 1 ) ±,,(A ^ A) , 

and 

MA)MA) = (A - A,)£(A - A s ) w (exp2i 7 (A)) , 
showing integrability of the singularity at A s . 

2. Estimation of P{X) and I q (X) 

In order to estimate / P (A) and I q (X) in the large-i limit, we use the steepest-descent method to evaluate the 
integrals (fB"5)) and ([B"i]) . The paths of the steepest descent going through the stationary point A s are shown in Fig.Q] 
An important consideration is that besides the contribution to the integrals of this stationary point, which is of the 
order 



one also has to take into account the contribution of the residues located at A ± iO and at the zeros of the function 
e A -/ 3 _ e ITrK . We begin by first neglecting the contributions from the residues at the zeros of e A _/3 — e' i7rK which at large 
t give exponentially small corrections and focus on the residue at A±iO. (A more complete estimate will be presented 
in the following sections.) Transforming the integration contour in (fSU)) from the real axis to the steepest-descent 
path T p (see Fig. [TJ, and using the analytical properties of the integrands discussed above, we obtain: 



7* (A) = V (X S A)MA)MA)p(A) e - 2 ^^ + O ( ^T^r ) 



(57) 
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FIG. 1: Stationary-phase contours for evaluation of the integrals (|53[) and (|54[) in the large-f limit in the space-like case. The 
dots are the zeros of the function e x _/3 — e 17VK with \q denoting the zero which gives the exponentially decreasing correction 
with the slowest rate of decay. 



Similarly, 

/£(A) = -„(A - \ s )5 + (\)S_(\)p(\)e-^^ + O ( —J— ) . (58) 

For I 9 (X), the computations follow the same steps with the steepest-descent path T q (see Fig. [1]), and the result is: 

4(A) = ±rj( T X s ± X)6^(X)5Z 1 (X)q(X)e^ t ^ + O ( —J-— ) . (59) 

The calculations above are valid when A is not too close to the stationary point A s . In the vicinity of the stationary 
point, the integrals I P (X) and I q (X) can be estimated as 



and 



with 



2m 



+°° (//-A.)yfc-A fl )irexp(-2iMQ) 
-oo /i — A 

+°° ( M - A,);"> - X S )Z W exp(2z0( M )) 

oo (J, — X 



dfj,, 



d/J., 



p s = p(X s ) exp(2ij(X s )) , q s = q(X s ) cxp(-2i7(A s )) . 

This means that in the vicinity of X s , the i-dependence of I p and I q in the leading order is given by the following 
relations: 



and 



+°° l M l 2 ^exp(-2it M 2 ) ^ = t _ iv f +°° l/x| 2 ^exp(-2i M 2 ) ^ ^ 

-oo A 4 - (^ - ^s) J-oo 



/i - \/i(A - A s ) 



+°° l A i|- 2ll 'exp(2zi/i 2 ) J _ +iv f + °° |/i|- 2w exp(2i M 2 



M-(A-A S ) 



-d/i = V 



(i - V~t(X - A s ) 



-dfi . 



(60) 



(61) 



The boundary values of these Cauchy integrals are uniformly bounded (see [38j) in i/i(A — A s ) due to the fact that v 
is real. This proves that the boundary values of I p and I q and therefore <&™(A) are bounded in the large-t limit. 
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B. Approximate solution of the RHP 

Now we are ready to show that the Manakov ansatz <& m (A) is an approximate solution of the RHP (|4"Tj) . More 
precisely, if ^(A) is the exact solution of (HJ), then 

<D(A) = [/ + 0(r^p m (A), for t^+oo, -±<-y/p, (3 > . (62) 

Indeed, not too close to A s , we have from Eqs. ([57]) , |58|) , and ([59)) : 



and 



§ »( A )=( i _ 1 ?7 (A-A s ) ( 5 + (A)A_(AMA)e- 2 ^) \ CT3 i n5 _(A) +0 



»7(A a -A)5^ 1 (A)*Z 1 (A)g(A) e 24 ' i W 1 J" \t/1(X — A s ) / ' 

so that 

,-2itA 2 -2i:rA 



[WCA)]" 1 *^) = ( .... '•.,„. p( , A)e '_ )+0 



g ( A ) e 2iU +2ixA l +p (A)g(A) J \y/i(X-X s ), 

Since the boundary values of (|60[) and (|61[) are uniformly bounded, this means that 

[$-(A)]- 1 $ m (A) = G,(A) + r (A) , (63) 

where G$(A) is given by (JUJ) and 

01 j tO(l), |A-A,| <t- 1 /2+e, ^ 
for t — > oo and with (0, -g)- If one introduces the matrix 

i?(A)^$(A)[$'"(A)]- 1 , 

then 

R+(X) - R_(\) = *+(*)[*+ (A)]" 1 - $_(A)[$ m (A)]- 1 , 

and using Eq. ([6"5[) and the relation $ + (A)G$(A) = $-(A) we obtain 

i? + (A) - R_(X) = r(A) = $ + (A)r (A)[$'"(A)]- 1 . 

Taking into account that i?(oo) = / we see that this this relation implies that the the matrix i?(A) can be represented 
like this 



1 rful 
RW = / + — / -^-du , for A e C/I 
2tt« J-oo M - A 



Under the hypothesis that $+(A) is uniformly bounded in A 6 K (which can be proved as in [33]), r(A) satisfies the 
same estimates as tq(A). Therefore, outside of a vicinity of A s , 



R(X) = l + 0(t-e), e e(o,i) 



proving Eq. 
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C. Asymptotic behavior of the potentials 



Making use of the Manakov ansatz 

$mm / -5(A) -I»(\)5- l M 

9 ~ y -I q {\)6{\) S-^X) 

one can extract the auxiliary potentials from the large-A expansion using the formulae obtained in Section IIV CI 
We start with b ++ which enters directly the expression (fTi?)) for the field correlator. Since <£> m (A) is an approximate 
solution of the RHP, Eq. can be written as 

1 r + °° 

b ++ = ($Di2 + o(l) = — J S+iriS^pirie-^^dn + o(l) , (65) 

with p(X) given by Eq. (|43|) . The integral appearing in this expression can be estimated via the steepest-descent 
method in the same way as we did for I P (X) in Section IV A 21 This means that if one neglects the exponentially small 
corrections that come from the residues at the zeros of e A _/3 — e 47rK , Eq. ([65]) gives 

b ++ =c a t- 1 / 2 -™e 2 ^ +o(l), 

where X s = —x/2t, v is defined by Eq. (|55p. and cq is a constant which depends on (3 and k. Until now, all the 
considerations were rigorous. The fact that the Manakov ansatz is only an approximate solution of our RHP, as 
specified by Eq. (|62p. means then that the next term in the asymptotic expansion could be of the order of 0(i _1 / 2_e ), 
and one should not take into account the exponentially small terms which appear in the complete evaluation of the 
integral (f65|) . There is, however, a caveat. The condition ([52]) ensures that the transformation of the integration 
contour from the real axis to the steepest-descent path encloses the pole at A = Aq which is closest to the real axis 
(see Fig. [1]) and is given by: 

A =-(/?+ VP 2 + t 2 k 2 ) V2 / V2 - i (-0 + v//? 2 + tt 2 k 2 ) 1/2 fy/2 . 

The residue at A = Aq~ gives the most slowly-decaying exponential term to the integral (|65p . i.e., a more complete 
solution for (|65|) including the contribution from A = Aq is 

b ++ = c t- 1/2 - 4!/ e 2ltA " +Cie - 2 ^(^,V) + (l), 

= cot-^ 2 - iv e 2itX ' + Cie 2*(™-;/3) e -2i*A - + _ (66) 

As one approaches the bosonic limit k = 0, the second term in (|66[) which arises from the pole at A = Ag~ becomes 
dominant even compared with a possible 0(t~ e ) term, since SAg~ — > in this limit. This term is the main component 
of 6++ in the case of impenetrable bosons. This shows that the exact solution for 6 ++ should be written as 

b ++ = c Q t- 1/2 - w e 2lt ^ +■■■ + Cie 2t(™-i0) e -2ix\- + . . . _ ( 67 ) 

where the dots between Co and c\ mean that there might be terms of order 0(t _1 / 2 ~ e ) which are, however, smaller 
than the c\ term when k — > 0. 

Although we will not use it below, we present the result for the potential B which is 

1 r +oa 

= + (i) = -—. J^M^'WpM^^'''"'^ + o(i) . 

This means that 

= c t- 1/2+lv e- 2ltx2 ° + o(l) . (68) 

By contrast, the results for potentials B^ and C_| + C |_ + B G will be very important for the subsequent 

calculations. They are: 

B+- = -oo - ($Dn + o(l) = -a - So + o(l) , 
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and 

C+- + G_+ + B__G = ($™) 22 - ($ 2 n ) n + o(l) = -25 x + o(l) , 
where ao is denned by Eq. (|48j) . and 

S = - dfi\n\ip(fi 2 ,f},K)\ , <5i = - / cfyt /xln |(^(/x 2 ,/3, k)| . 

7T / ™ 7T /_„„ 



Introducing the function J(/3, k): 

r + oc 



I(P,K)= $slj d[i In ip(p 2 ,/3,K)j , 

we can rewrite the result for B_| as 

B+- = + — / sign( M - A.) In |^(/i 2 , 0, k)| dp + o(l) . (69) 

Zir 2tt J_ oq 

Also, using the fact that <p(/i 2 ,P, k) is even function of p, one can transform the result for C-\ + C |_ + 5__G into 

A r + oc 

C+- +C-+ + G = - / sign(p- A s )pln|<p(p 2 ,/3,«0| dp + o(l). (70) 

D. Asymptotic behavior of a(x,t, f3, k) 

The formulae ()69|1 and (I70|) allow us to obtain the asymptotic expression for a. As a first step, combining them 
with the differential equations (f2"6")) we have: 

x 1 



r+oo 

d x <r=-I(fi,K) + - I sign(p- A. s )ln|^(p 2 ,/?, K )| dp + o(l) , (71) 

7T 7T 



and 



2 



+ 00 



d t <T=- sign(fj,-X s )fj,hx\(p{fJ,\0,K)\d f i + o(l). (72) 

The asymptotic expression for a is obtained integrating Eqs. (fTTj) and (f?2")) over a; and i. This implies, however, that 
more accurate expressions for the derivatives of a that include the higher-order asymptotic terms are needed to have 
the same accuracy for a as for b ++ (|6"T|) . To obtain these expressions we first note that Eq. implies that 

d x B + _ = ^ln\<p(\ 2 s ,j3,K)\j = o(^) . (73) 

Combined with the first part of Eq. (|24| . d x B^ = 2i6 ++ _B__, this result agrees with the estimates b++ = Oil/t 1 / 2 ) 

and = 0(l/tV 2 ) that were already obtained in the previous section, see Eqs. (|66|) and (|68|) . Also, we know that 
the potentials and b ++ solve the separated nonlinear Schrodinger equation (|2"3")l for which the general structure 
of the decreasing solutions is (see, e.g., [H, 

u a + 2^1^ — 7^ — I e , (74) 

n=l fe=0 / 

B- = t-" 2 ( ,0 + E E ^5^1 e^^ 1 ^ , (75) 
\ n=l fc=o / 
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where uo,vo,u„k,v n k and v are functions of A s = —x/2t. The parameters v,u n k,v n k can be expressed in terms of 
uq, Vq. In particular, 



Vi 2 Uq + U12W0 = , v = ~Au v , 
(„»)» 



v w u + u 10 v 



32 ' 

{vu'Y i , ... 



(76) 



where the prime denotes the derivative with respect to A s . Now we can improve the asymptotic expansions for the 
derivatives d x a and dto. Substitution of (f74|) and (j75|) into d x B^ = 2ib ++ B gives 



d x B, 



iv (v 2 )"\nAt 



2\ll 



t 2 \ t 3 



(77) 



Comparing the first term in this expansion with Eq. (|T3"|) . we see that v = — — In \f(X 2 , 0, k)\ > in Eqs. ([74")) and 
(I75p . in agreement with our previous notation (|55|) . Integrating Eq. (|77p over cc and using the first equation in (|26[) . 
we obtain 



j 1 f +oa 

d x a = -1(0, k) + - / sign(fi- \ s )\n\<p(n 2 ,0,K)\ dfi 
n ^ J-oc 



(^ 2 )'ln4i 



(v 2 )' 



+ i(v' u Q - vqUq) 



4 t 

Equation ([23]) can be rewritten as 

d x (C+^ + C-+ + B__G) = b ++ d x B__ - B__d x b ++ . 
Then the asymptotic expansions (|74[) and ([75)1 give 

^^((j/ 2 )') xd x (v' u - u' v Q ) (u v' - v u' ) 



d x (C + ^+C^+ + B^G)=i- 2t2 



sr 2 



2f- 



2t 2 



(78) 



xdMv 2 )') ln4i ^(9 x ^)ln4t ^/ln 4 4i 
+ i — — — — — - i K ; ' — — + O 



8 t 2 2 t \ P 

Integrating this equation over x, and using the second equation in (|2"6")l we find 



d t a = - 

IT 



sign(/z- A s )/^ln |<^ 2 ,/3, k)| dfj,-\ 6 



2 i + 2* 



(79) 



A s -^ X s j(v u - v u a ) + O I — — 



(80) 



Finally, integration of Eq. (|79[) over a: and (|B0[) over t gives the asymptotic expansion for a(x,t, 0, k) of required 
accuracy: 

i 1 /"+ 00 i/ 2 j/ 2 

a(x,t,0,n) = x-I(0,k) + - / \x + 2t\\n\p(/j 2 ,0, k)\ -\ 1 In At 

7T 7T J_ OQ 2 2 



2i / ( v 'o(n)uo(n) ~ v (n)u' (p) dfi + c(0) + O 



In 4 At 



(81) 



where c(0) is a constant that depends only on 0. 
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VI. ASYMPTOTIC SOLUTION OF THE RHP. TIME-LIKE CASE 

The computations in the time-like region defined by 

\ S >~^P, = h/T>O, 

are very similar to those presented above for the space-like case. Because of this, the presentation in this Section is 
more sketchy, emphasizing the differences between the two regions. As we will see in what follows, the leading term 
of the asymptotics for the potential 6 ++ is the same in the time-like as in the space-like region. The sub-leading term 
in the asymptotic expansion, which is important because it reproduces the predictions of conformal field theory, is, 
however, different in the time-like case. 

A. Manakov ansatz 

The Manakov ansatz in the time-like region is: 

4> m (A) = ( _ J J (A) \ e -3in*(A) ? (82) 

where / P (A) and I q (X) are now given by 

P(X) = -L / + °° ^'^ -^ , (83) 
2tt« J-oo M ~ A 

and 

i«(A) = -1- [ + °° t+W-yi qUyM*'**)^ . (84) 
2tt« J~oo M - A 

The functions p(A) and q(A) here are defined in Eqs. (|46|) and (|47|) . respectively, while the function 8(A) is the solution 
of the following scalar RHP 

6+{X) = 6-(\)[l+p(X)q(X)r}(\-\ s )], Ael, 5(oo) = l. 

Using the fact that for p(X) and q(X) defined by (l4"6")) and l|4"7|). l+p(X)q(X) = \ip(X 2 ,(3, — k)\ 2 , one can see that solution 
of this RHP can be written as 

*«=e*p{^/J^l» W( AA-«)l>}. 

1. Properties of 5(A) 

Following the same steps as in the space-like region, we have 

MA) = (A - W exp(i 7 (A))KA2,/3, -«)|± (|^(A 2 , 0, - K )|) ±7)(A - Xs) , (85) 

with 

l(A s ,/3, k) = ~]n\<p(\ 2 a ,p, k)\ = --In \p(A 2 s , (3, -k)\ > , 

7T 7T 

and 

7 (A) = - / ln|/i-A|d(ln|^( M 2 ,/3,-«)|)d/x. 
A, 

This also means that 

S+(A)S.(A) = (A - A S )+* V (A - A s )Z w (exp2z 7 (A)) , (86) 
showing integrability of the singularity at A s . 



19 




FIG. 2: Stationary-phase contours for the integrals (|83|l and (|84l) in the large-t limit in the time-like case. The dots are the 
zeros of the function e A _/3 — e~' 7VK with Aq~ denoting the zero which gives the exponentially decreasing correction with the 
slowest rate of decay for &++ in the time-like case. 

2. Estimation of P{\) and I q (X) 

The estimates of I P (X) and I q (X) in the time-like region are obtained as in the space-like region by the steepest- 
decent method. The steepest-descent contours are shown in Fig. [2j Similarly to the space-like case, for A not too 
close to the stationary point A s , transformation of the integration contour from the real axis to the steepest-descent 
paths gives: 

4(A) = ± V (±\ S T X)5^(X)SZ 1 (X)p(X)e- 2i '"^ + O ( - , (87) 



x s )J' 



and 



4(A) = ±r,( T A s ± X)5 + (X)5-(X) q (X)e 2 ^-^ + o ( — J— - ) . 

Also similarly to the space-like case, one can show that for A in the vicinity of the stationary point, the boundary 
values of I p and I q , and therefore $™(A), are bounded in the large-t limit. 

B. Approximate solution of the RHP 

Combining Eqs. j87j), (88]), and (82]), one obtains 

ram/^-i _^m« + (A) ( 1 r 1 {X s -X)5 + {X)-^Z 1 {X)p{X)e-^W \ ( 

L +WJ ~ & { V (X-X s )5 + (X)5 ( X)q(X)e^W 1 ) + ° \ 



and 



*m m / 1 7 7 (A-A s ) ( 5 + (A)- 1 A: 1 (AMA)e- 2 ^ A )\ CT3ln5 _ (A) / 

- {X) -{v(X s -X)S + (X)S^X)q(X)e^W i ) e +°{ 



Vt(x - x s ) 



\/i(X-X s ) 



and therefore, 



[$+(A)] $ - (A) -U(A)e 2 ^ 2+2 ^ 1 ) + °{Vt(X-X s ) 
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This shows that the Manakov ansatz ([82)) is an approximate solution for the RHP (|4Tj) with the conjugation matrix 
P5| . More precisely, if ^(A) is the exact solution, then 

$(A) = [J + 0(t^)]$ ro (A), eefaiY for i^+oo, -| >-v^, /? > . 



C. Asymptotic behavior of the potentials 

As above, the asymptotic expressions for the potentials are extracted from the large-A expansion of Eq. (|82|) making 
use of the formulae obtained in Section HV CI Substituting Eqs. (|82|) and (|83|) into the second equation in (|49|) . we 
have 

b ++ = ($r)i 2 + = t— / 5 + ( M )- 1 5-(M)- 1 p(M)e- 2#(x '^ ) ^ + , 

where p(X) is given by Eq. (|46p . The main difference with the space- like region is that now the residues that give 
the exponential corrections to the leading term of the asymptotics are the zeros of the function e A ~^ — e~ l7rK and 
not e A ~P — e l ^ K . The pole that is closest to the real axis among those that are enclosed by the transformation 
of the integration contour from the real axis to T p (see Fig. [2]), and therefore contributes the most slowly-decaying 
exponential term, is: 



A, 



(B + V/3 2 +7T 2 K 2 ) V2 /V2 + i C-j3 + V/? 2 + 7T 2 K 2 ) V2 /V2 



(I 



The contributions of the stationary point and the residue at Xq produce then the following expression for b ++ : 

b ++ = c t- 1/2 - w e 2ltx2 s + Cie -2^(x,j/,A -) + ( 

= C t" 1/2 ~ W e 2itA ' + Cie 2t(-T«-i« e -2foA- + ^ (gg) 

where Co and ci are some undetermined amplitudes which can depend on /3, K, and A s . Again, as we approach the 
bosonic limit, k — > 0, the second term in (|89[) becomes dominant. This term represents the leading asymtoptic term 
of 6 ++ for impenetrable bosons. 

The leading term of the potential B is given by the same Eq. ([68]) as in the space-like case. The potentials 

and C H + C |_ + B G are obtained from Eq. ([50|l 

B+- = +a - ($r)n + o(l) - +a -S + o(l) , 

and 

C+- + C-+ + B__G = ($5% - ($ 2 n ) n + o(l) = -25 x + o(l) , 
where o/q is now defined by Eq. (|5ip . and 



If 00 1 
So = — dfiln |tp(/i 2 , (3, -k)\ , Si = — I dfJ, /Ltln |<p(/i 2 , ^, -k)| . (90) 



7U 

The result for Bj can be rewritten as 



B+- = J(/3 ' K) + ^ / + °° sign(/i - A s ) In |^ 2 , /?, -«)| d M + o(l) . (91) 



>+- 

2vr 2tt _ 

in terms of the function 

I(/3,—K) = ^i / d/i In </?(/j, 2 , /3, — k) I = — Q I / d/i lrnp(fi 2 , f3, k) 

Also, using the fact that (p(/j, 2 , f3, —k) is an even function of n, the equation for C_| + C |_ +B G can be transformed 

into 

■ p + OC 

C+_ +C_+ + £__G= - / sign(Ai- A s )^ln |^( M 2 , /3, -k)| d/x + o(l) . (92) 
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D. Asymptotic behavior of a(x,t, /3, n) 

Taking into account that 1(0, —k) = —I(f3, k), and \tp(n 2 , 0, —k)\ = |<p(/i 2 , 0, one sees directly that the asymp- 
totic expressions (|9"T)l for i? H and (152")) for C + _ + C |_ + £?__G in the time-like case coincide with the corresponding 

expressions and ([70)1 in the space-like region. Since the higher-order corrections discussed in Sec. IV Dl are the 
same in both regions, this means that the asymptotic expansion for a in the time-like case is given by the same 
Eq. (JHU) as before. 



VII. RESULTS 



Now we have all the ingredients to formulate the results for the main object of our interest, the anyonic field 
correlator which, as a reminder, is given in rcscaled variables by the expression 

(*(as 2 , *2)* t (zi, t x )) T = VTg(x, t, 0, k) (93) 

with 

g(x, t, 0, k) = -^-e 2it ^b ++ (x, t, (3, K ) e ^M . (94) 



A. Negative Chemical Potential 



While all the considerations in the previous sections were based on the assumption that the chemical potential is 
positive, in fact, the results obtained are also valid when (3 < 0. In this case, we need only the leading term for b ++ . 
Putting together the first term in Eq. (JSTJ or ([55| and Eq. (I5T1) . we can express the leading asymptotic behavior of 
the anyonic field correlator at negative chemical potential as 




where cq is some constant amplitude, A s = —x/2t, v = —(1/tt) In \tp(X 2 ,0, k)\, and the definitions of all other functions 
in this equation are presented together in Eqs. (J5J) and © of the Introduction. 



B. Positive Chemical Potential 



For reasons discussed in Sec. IV CI in the case of positive chemical potential, one needs to keep in the asymptotic 
expansion of the potential 6++ and, therefore, of the field correlator, not only the leading term, which is the same in 
the space-like and time-like regions, but the next exponentially decreasing term as well, which is different in the two 
regions. Thus, the two results should be presented separately. 

Space-Like Region: x/2t > \/]3. Combining Eqs. (|8T)) and ([67| . we have 



g(x,t,0,K) = ^ 2 /2 e «/(/3,K)A e C(z,M,K) ^ CQt -l/2-iv e 2it(\l+0) +Cie 2t™ e -2i X X; +Q f t -l/2 

Time-Like Region: x/2t < \/]3. In this case, Eqs. (|8T)) and ([89| give: 

g(x,t,/3,K) =t « 2 /2 e i*WK)/* e C(x,t,P,K) ^ Cot -X/2-iu e 2it(Xl+P) + Cie -2t™ e -2«A+ + Q ^-1/1 

The constants and Xq in these equations are defined by Eq. §5§ of the Introduction. 



(96) 



(97) 



VIII. BOSONIC AND FREE-FERMIONIC LIMIT 



As the last step, we analyze our main result for the anyonic field correlator in various limits, in order to establish 
the relation with previously known expressions, and to demonstrate the unexpected features of the anyonic case. 
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Bosonic Limit. For bosons, k — > 0, one has: 

A 2 -/5 _ i i H-oo 

<p(\ 2 ,f3,K = 0) = -r^— = , C(x,t,0,K = 0) = - / |x- 2iA|ln|(^(A 2 ,/?,K = 0)| dA, (98) 

e A p + 1 ttJ-oo 

and i/ = - (1/tt) In | <p(\ 2 s ,0,K = 0)|. Also, 1(0, K = 0) = -2ttV? for /? > 0, and 1(0, K = 0) = for < 0. In the case 
of negative chemical potential, using these relations it is straightforward to see that Eq. ((95]) reduces to the known 
result for impenetrable bosons 0, [33| • For positive chemical potential, the result obtained in [H, [33| is 

g(x, t,0,K = O) = Co t v2 l 2 e c{x ^' K =^ \l + 0(t-^ 2 )] , (99) 

and is valid in both the space-like and the time-like region. Taking into account that in both regions, A^ = —\[0~ 
for k — 0, we can see that in Eqs. (|96|) and (l97|) . the second term in the parenthesis gives the leading contribution 
in this limit, which reproduces the bosonic result. This means that for a certain value of k approaching 0, there is a 
crossover in which the relative magnitude of the two terms in the parenthesis changes, and the second term becomes 
the leading one for k close to 0. 

Free-Fermionic Limit. For k — > 1, the anyonic system we considered reduces to free fermions. In this case, the 
function (p(\ 2 ,0,K) vanishes, which means that v = and C(x,t, 0, k = 1) = 0. It is easy to see that Eqs. (|96p and 
(|97|) reduce to the corresponding correlators (|A1[) and (|A2[) of free fermions that are presented in Appendix fA"l 



IX. CONFORMAL FIELD THEORY 



The behavior of the field-field correlators of the one-dimensional particle systems at low temperatures is usually 
believed to follow the predictions of conformal field theory (CFT). For impenetrable anyons, the CFT result for the 
leading term of the large time and distance asymptotic of the field- field correlator is [3] : 

mix, mU0, 0)) - e~ ikFKX exp (- ( ^-^—\x - v F t\ + ^^—\x + v F t\ 

I \ V F V F 

where k F = \fh and v F = 2k F are the Fermi vector and the Fermi velocity, respectively, and A ± are the conformal 
dimensions: 

(Note the differences in conventions between Q and this work together with the rest of the papers [l|, [|J , where we 
have obtained the determinant representation for the anyonic field-field correlator. The main difference is related to 
the ordering of the anyonic creation operators in the eigenstates of the Hamiltonian, which amounts with the change 
of the sign of the statistics parameter k in the formulae of Q.) In the notations of this work, the field correlator 
considered here is 

{^(x 2 ,t 2 )^(x 1 ,t 1 )} T , t 21 =h-t 1 >0, x 12 =x 1 -x 2 >0 } 
which means that the CFT predictions are 

{^(x 2 ,t 2 )^(x 1 ,t 1 )) T - e ikFKXl2 e nTKt2i e~^ Xl2 ( i+ ^), (100) 

in the space-like region, and 




(*(ara,*2)* t (a;i,ti)>T ~ e ^ Ml2 e^ 12 e"" Tt21 ( l+ ^), (101) 

in the time-like region. The leading term of the asymptotics (|96|) and (|97|) obtained from the exact calculation in this 
work do not reproduce these equations in the limit of low temperatures — * oo. We can show, however, that the 
sub-leading terms in these asymptotics do give the conformal behavior. Indeed, as shown in the Appendix [B] in the 
limit — > oo, we have: 

x-I(0,k) +C(x,t,0,K) = x2iJ0U - l)-x^-=(l - k) 2 , (102) 

7T 2 VP 
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in the space-like case, and 

x-I(fl, k) + C(x, t, (3, k) = x2i^/f3(n - 1) - *7r(l - k) 2 , (103) 

7T 

in the time-like case. In the same limit, Xq and Xq given by ([5]) become 

X± = -^3 ±l ^=. (104) 

Using these formulae, we see directly that the second term in the asymptotic expansion of the field correlator is given 
by 

in the space-like, and 

r 2 /2 e 2 l ^v? e ^ e -t-(« 2 +i) j (106) 



in the time-like region. The exponential terms are exactly the ones predicted by CFT, if we take into account that 
x = xisVT, t = t 2 iT/2, and /3 = h/T. 

Qualitatively, the non-conformal term of the time-dependent field-field correlator, which is the leading asymptotic 
term for particle statistics not too close to bosons, can be traced back [4l[ to the singularity of the one-dimensional 
density of states at the bottom of the single-particle energy spectrum A — > 0. In agreement with this interpretation, 
there is no non-conformal terms in the "static" equal-time correlator (see, e.g., [29|). since the single-particle spectrum 
is unlimited in momentum space, A E (— oo, +oo). By contrast, the energy spectrum e oc A 2 has a threshold at A = 
with associated non-analytical behavior of the density of states. This non-analyticity manifests itself directly through 
the non-conformal terms in the asymptotic behavior of the field correlator of the massive one-dimensional particles. 



X. SUMMARY 



In conclusion, we have calculated the large time and distance asymptotic behavior of the temperature dependent 
field-field correlation functions of impenetrable one-dimensional spinless anyons. As a function of the statistics pa- 
rameter, the anyonic correlator interpolates continuously between the two limits of impenetrable bosons and free 
fermions. The main qualitative feature of our result is that, asymptotically, the anyonic correlator consists of two 
additive parts. One is a non-conformal term produced by the non-analyticity of the density of states at the bottom 
of the single-particle energy spectrum. For all values of the particle statistics away from the bosonic limit, this term 
gives the leading asymptotic contribution to the correlator. The other is the sub-leading term which agrees with the 
conformal field theory and is associated physically with the low-energy excitations close to the effective Fermi energy 
of the system of impenetrable anyons. In agreement with the previous results 0, [4^], for the statistics parameter 
close to the bosonic limit, the conformal term determines the leading behavior of the asymptotics. Because of the 
additivity of the two parts of the correlator and their different physical origin, even away from the boson limit, the 
systems response to the low-energy probes is determined by the (sub-leading) conformal part of the correlator. 



APPENDIX A: LARGE TIME AND DISTANCE ASYMPTOTIC BEHAVIOR OF THE FIELD 

CORRELATOR FOR FREE FERMIONS 



In rescaled variables used in this work, t = {t% — t\)T/2 > 0, x = (x\ — x<x)\[T /I > 0, (3 — h/T, the field-field 
correlation function of free one-dimensional fermions is expressed as [42[ 

2itp r+oo p \ 2 -(i 

(*(z 2) t 2 )^i,ii)>T = VT— J dX eX2 _ p + i e^f^\ 

We are interested in the asymptotic behavior of the correlator in the limit of large x > 0, t > with x/t — const. 
The analysis is similar to the one performed for the functions / P (A) in Section TV A 21 The leading term is obtained 
via the steepest descent method and the corrections come from the poles located in the complex plane at the zeroes 
of the function e A + 1. The corrections to the leading term are different in the space- like and time- like regions. 
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Space-like region: (x/2t > \f]5). In this case, the residue that gives the exponential term with the slowest rate 
of decay for large x and t is 

X s = - (p + ^WT^ 2 ) V2 /V2 -i (-13+ ^fWT^) ^ /V2 , 
resulting in the following asymptotic behavior of the correlator: 

(*(x 2 , t 2 )^(x 1 , h)) T ~ cot- 1/2 e 2ii ^+ A ^ + Cl e 2 * T e - 2 " A o + • • • . (Al) 

Here A s = —x/2t is the stationary point of the phase <f>(x,t, A) and cq, c\ are some constant amplitudes. In the limit 
of low temperatures, (3 — *■ oo, using the fact that Aq — ► — V/3 — «7r/ (2y/]3~), we obtain 

(tt(z2,ta)* + (a:i,ti)>T ~ c t- 1 / 2 e 2 "^+ A ") + Cie 2 ^ e M*-/ 2 v?) + . . . . 
Time-like region: (a;/2t < \f]3). In this case, the residue producing the leading contribution is: 

\l = - (j3 + ^WT^ 2 ) 1/2 / V2 + 1 (-p + ^WT^ 2 ) 1/2 / V2 . 

with the corresponding asymptotic behavior of the correlator: 

(^(x 2 ,t 2 )¥(x 1 ,t 1 )) T ~ cot^e 2 ^^ + Cl e- 2t *e- 2txX ° + ■■■ . (A2) 
In the low-temperature limit, this becomes 

(^(x 2 ,t 2 )^(x 1: t 1 )) T ~ c t- 1 / 2 e 2W ^+ A ^ + Cie 2 ^^ e 2 -(-*+-/ 2 ^ + ... . 

APPENDIX B: ANALYSIS OF C(x,t,f3,n) 

The function C(x, t, (3, k) is defined in the main text as 

1 f +DO 

C(x,t,/3,K) = - \x-2t\\\n\<p{\ 2 ,(3,K)\ dX, (Bl) 

7T ./_„ 



where 



y(A 2 ,/3, K ) = + l • (B2) 

Using the expansion of the logarithm: ln(l — z) = — X^^Li z n /n , |z| < 1, we obtain the following expansions for 
ln|^(A 2 ,/3, K )|: 



ln|^(A 2 ,/3, K )| =-J2 (cos(mm) + , Ae(-^,^), (B3) 



e 

n=l 

and 



°° p n(/3-A 2 ) 

ln|^(A 2 ,/3,«)| =- V (cos(n7r«) + , A G (— oo, —y/P) U °°) • (B4) 

n=l 

We are interested in the the asymptotic behavior of C(x,t, (3, k) in the limit of low temperatures (/3 — * oo). This 
behavior is different in the space-like and time-like region. 

Space-like region: (x/2t > \f]3). It is convenient to express the function C(x,t, j3, k) in this case as 

i rx/ 2 t 1 />oo 

C(x,t,(3, K ) = - (x-2tX)ln\(p(X 2 ,(3,K)\ dX + - {2tX - x) In |<^(A 2 , 0, n)\ dX . (B5) 

7T J-oo n Jx/2t 

Using the expansion (|B4[) one can see that the second integral in this equation is on the order of O ^e - ^/ 2 *) 2- ^ , 
which for x/2t outside of the immediate vicinity of y/]3, more precisely: x/2t—\f]3 > 0{l/\fjJ), decreases exponentially 
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in v 7 /?, since (x/2t) 2 — (3 > 2y/j3(x/2t — \/~j3). The same argument allows us to extend the upper limit of integration 
in the the first integral on the RHS of Eq. (|B5[) back to +00. Then, the expansions (|B3[) and (|B4[) combined with the 
formulae 

e /" e x2n dX = +o( , [°° e~ x * n dX = — ^ + O ( ' 



2n^]3 V/? 3/2 /' -A/3 2nV^ V/? 3/2 /' 

give the following estimate for this integral: 

cos(roivc) + 



i / (x - 2tA) In \ V {X\ 0, K )\dX = -x* T, C ° S(n7rK) t + O 

T '-co ' "'VP ^ ™ 



£3/2 



(B6) 



Using the formulae (0.234) and (1.443) of |43|: ££l 1 (-l) n+1 / n2 = ^V 12 and J2T=i cosnirn/n 2 = tt 2 B 2 (k/2), where 
-62(2;) = x 2 — x + l/6is the second Bernoulli polynomial, and the fact that contribution of the region A > x/2t to the 
integral can be neglected, we rewrite the previous result as 



\ x -2tX\ln\<p(X z ,P,K)\dX = -x^(l- K y+0\-^), (/3 - 00) . (137.) 



Therefore, in the space-like region, we have 
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C(x,t,p, K ) = -x^=(l-K) 2 +Ol^^y (0 -00). (B8) 

Time-like region: (x/2i < In this case, we begin by expressing C(x, t, f3, k) as 

4 r°° 2 
C(x,t,{3,K)=t~ Aln|v?(A 2 ,/3,K)| dX- - / (2At - x) In |^(A 2 , 0, k)| dX . (B9) 



As above, using the expansion (|B3[) . one can see that the second integral in this equation is on the order of 
0(e~(P~( x ' 2t ' which for x/2t not too close to \/]3, more precisely: x/~j3 — x/2t > 0(1/^/0), decreases exponen- 
tially in y/]3, since (3 — (x/2t) 2 > 2\//3(\/l3 — x/2t). The, the expansions (|B3|) and (|B4|1 and the calculations similar 
to those in the space-like region give for the first integral: 

4 f°° 

t- Xln \<p{X 2 ,f3, K )\ dX = -irt(l - k) 2 , (BIO) 

7T Jo 

The final result is 



C(x, t, (3, k) = -7rf(l -nf + O (e-W-W'A • (Bll) 
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